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Abstract 


Snow studies drastically influences on managing water resources and preventing the dangers caused by floods. 
Moreover, the modeling of snow melt runoff is faced with the lack of snow density data in the basin due to many 
changes in the distribution of snow masses, and creating snow measurement stations at high altitudes are 
difficult and expensive work. Considering the issues, this study attempted to invent a new equation and use 
remote sensing techniques in a simpler and more physical way to determine the degree-day factor (a) and then 
the simulation was performed by defining its calculated values using classical and new methods along with the 
annual and seasonal values of the recession coefficient (k) for the snow-melt runoff model (SRM). The 
outcomes were evaluated and compared. In order to evaluate the performance of the model and its parameters, 
simulations were performed for the calibration and validation periods for the water years 2011-2012 and 2012- 
2013, respectively, and the MODIS snow cover product (MOD10A1) was used to estimate the snow-covered 
area and to estimate the radiation values. NEO net radiation product (CERES-NETFLUX-E) was used. The 
results showed that the best method for simulating runoff in the calibration period (2011-2012) was the new 
degree-day factor calculation method and the two values of the seasonal recession coefficient with the 
coefficient of determination of 0.72 and the volume difference of 4.17. In the validation period (2012-2013), the 
runoff simulation method with the new degree-day factor method and an annual recession coefficient with a 
coefficient of determination of 0.51 and the volume difference of 4.38 provided the best results in terms of 
assessment of the model accuracy criteria. 
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چکیده 

با توجه به تأثیر مطالعات برف در مدیریت منابع آب و جلوگیری از مخاطرات ناشی از 
بروز سیل و ازآنجاکه مدل‌سازی رواناب ذوب برف به دلیل تغییرات بسیار در 
پراکندگی توده‌های برف با کمبود یا فقدان داده چگالی برف در حوضه ppl‏ مواجه 
بوده و Bem‏ ایستگاه‌های اندازه گیری برف در ارتفاعات کاری سخت و پرهزینه است» 
در این مطالعه تلاش شده با ابداع معادله‌ای tte‏ و با استفاده از تکنیک‌های 
سنجش‌ازدور به روشی ساده‌تر و فیزیکی‌تر پارامتر فاکتور درجه-روز (0) را محاسبه 
کرده و سپس شبیه‌سازی با تعریف polis‏ محاسبه شده آن به روش‌های کلاسیک و 


جدید در polis LS‏ سالانه و فصلی ضریب فروکش K)‏ برای مدل رواناب ذوب 
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برف (SRM)‏ انجام شده و نتایج حاصل مورد ارزیابی و مقایسه قرار گرفته است. به- 

منظور ارزیابی عملکرد مدل و پارامترهای Ol‏ شبیه‌سازی برای دوره‌های واسنجی و 

درستی‌سنجی به‌ترتیب برای سال‌های آبی ۲۰۱۱_۲۰۱۲ و ۲۰۱۲۲۰۱۳ انجام شد و از 

محصول پوشش برف MODIS (MODIOAT)‏ جهت برآورد سطح پوشش برف و 

برای برآورد مقادیر تابش gale‏ از محصول تابش خالص NEO (CERES-‏ 

NETFLUX-E)‏ استفاده گردید. نتایج OUS‏ داد بهترین روش شبیه‌سازی رواناب در 

دوره واسنجی (۲۰۱۲-۲۰۱۱) استفاده از روش محاسبه فاکتور درجه-روز جدید و 

کاربرد دو مقدار ضریب فروکش جریان فصلی با ضریب تبیین ۰/۷۲ و درصد اختلاف 

حجمی ۶/۱۷ است. در دوره اعتبارسنجی (۲۰۱۳-۲۰۱۲) نیز شبیه‌سازی رواناب با 

روش محاسبه فاکتور درجه-روز جدید و یک ضریب فروکش سالانه با ضریب تبیین 

۱ و درصد اختلاف حجمی 1/۳۸ بهترین نتیجه را از نظر ارزیابی معیارهای دقت 

مدل ارائه کرد. 

کلیدواژه‌ها: ذوب برف. فاکتور درجه-روز» سنجش‌ازدور MODIS‏ تابش خالص 
۱- مقدمه 

ایران سرزمین کم‌آبی است که در منطقه خشک و نیمه‌عشک کره زمین قرار گرفته است و لذا بارش برف 
مهم‌ترین عامل تأثیرگذار در میزان ذخایر آب به‌ویژه در فصول گرم سال محسوب می‌گردد (قنبرپور و 
همکاران ۱۳۸۶). برف یک شکل از بارش است که به دلیل تأخیر در OF bas‏ به رواناب تفاوت ویژه‌ای با 
دیگر اجزاء دارد (عالی‌نژاد و همکاران» Y‏ رواناب حاصل از ذوب برف به دلیل نقش تأخیری آن منبع 
Lol‏ تغذیه سفره‌های آب زیرزمینی و در برخی از موارد به دلیل هم‌زمانی با بارش‌های oly‏ منشأً بروز 
سیلاب‌های مخرب با حجم جریان بیش از ظرفیت رودخانه‌ها می‌شود (میریعقوب‌زاده و قنبرپور» ۱۳۸۹). 
برآورد دقیق آب معادل برف و همچنین تغییرات سطح پوشش آن. در بخش‌های کشاورزی, انرژی؛ 

مدیریت مخزن و هشدار وقوع سیل مؤثر است (جهان‌بخش اصل و همکاران. ۱۳۹۵). شبیه‌سازی رواناب 
ذوب Gy‏ در مناطق کوهستانی نه‌تنها برای ارزیابی مخاطرات سیل ناشی از ذوب برف در بهار و QUA‏ 
بلکه برای توسعه و مدیریت ale‏ آب در حوضه اهمیت زیادی دارد (هاو و همکاران» ٩۲۰۱؛‏ اما در 
کشورهای درحال‌توسعه و مناطق دورافتاده با کوه‌های مرتفع» برآورد برف روی سطح زمین به دلیل 


y ۰ 2 -‏ 
توپوگرافی پیچیده عدم دسترسی و محدودیت‌های مالی بسیار سخت و محدود است (ویوویچ و جاکوبز » 


1 Hao et al. 
2 Vuyovich et al. 


سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل SRM‏ و برآورد پارامتر فاکتور درجه-روز با استفاده از... Y‏ 


QM‏ با وجود تغییرات زیاد در توزیع Sk;‏ و مکانی برف انباشته شده. اطلاعات مشاهده‌ای زمینی برف 
اندک و معمولا به نواحی با ارتفاع کمتر محدود می‌شود که این موضوع کاربرد سنجش‌ازدور در تأمین بخشی 
از اطلاعات برف را اجتناب‌ناپذیر نموده است. بعلاوه مدل‌های پیش‌بینی رواناب حاصل از ذوب برف نیاز به 
اطلاعاتی دارند که این اطلاعات با فقر آمار در مناطق برف‌گیر حوضه‌های آبریز عمدتاً قابل‌دسترس و یا 
تأمین نیستند که با این توصیف تکنیک سنجش ازدور به‌عنوان ابزاری جدید و مفید برای به‌دست آوردن داده- 
های برف مطرح می‌باشد (نجف‌زاده و همکاران. ۱۳۸۳). اندازه‌گیری‌های میدانی که به‌ویژه در مناطق 
کوهستانی انجام می‌شود اغلب بدون استفاده از فناوری‌های جدید غیرممکن است (شهابی و همکاران. 
۶ فناوری سنجش‌ازدور با آماده‌سازی لایه‌های اطلاعاتی پیوسته و با دقتی SVL‏ و هزینه کمتر می‌تواند 
Se‏ آمارهای هیدرولوژیکی را پر کند (طالبی اسفندرانی و همکاران. ۱۳۹۰). به‌طور کلی» مدل‌سازی ذوب 
برف و رواناب حاصل از Of‏ با استفاده از دو روش انجام می‌شود: روش تعادل انرژی و روش درجه-روز دما 
که روش دوم می‌تواند در GLE‏ داده‌های کافی مفیدتر باشد (سنزبا و همکاران؛ ۲۰۱۵). پارامترهای برف نقش 
بسیار مهمی در چرخه هیدرولوژیکی, الگوی آب‌وهوای محلی و تغییر اقلیم ایفا می‌کنند و درنهایت بر کل 
رواناب جریان رودخانه‌ها و آب در دسترس از طریق مقدار و زمان ذوب برف تأثیر می‌گذارند (نانچهانی و 
allis‏ ۲۰۲۰). 

در مرور پژوهش‌های مرتبط انجام شده در این adam‏ مارتینک (۱۹۸۹) در سایت آزمایشگاهی 
weissfluhjoch‏ واقع در آلپ سوئیس میزان ذوب ساعتی برف را در سراسر یک فصل ذوب برف محاسبه 
کرد و بده پوشش برف به‌طور مداوم توسط یک لایسیمتر برف اندازه‌گیری شد. palie‏ به طی تمام دوره 
فرسایش بیشتر در محدوده ۰/۲ تا ۰/۲۵ به‌دست آمد و مقادیر محاسبه شده می‌تواند برای مطالعات > CS‏ 
آب در پوشش برف استفاده شود. کوستاس و رانگو (۱۹۹۶) یک ad js‏ تابش ساده را با روش درجه-روز در 
ER T NE‏ بیرف یل IETNCETE FORM‏ مطالته: Seal‏ اکور ذوت 
برف فیزیکی‌تر برای SRM‏ توسط ترکیب روش بودجه تابش ساده شده با روش درجه-روز بوده است. برای 
روش درجه-روز محدود مقدار ثابت ۰/۲ برای » در طول فصل ذوب برف استفاده شد؛ در حالی که در 
روش درجه-روز bol‏ مقادیر » به‌تدریج از ۰/۶۸ در ماه مه افزایش CL‏ و به ۰/۵ در ژوئن و ۰/۵۲ در ماه 
ژوئیه رسید. بروبکر و همکاران" (۱۹۹7) پژوهشی در حوضه آبریز اسلیپرز ریور انجام دادند. هدف از این 


مطالعه ol ($2555 —]Á‏ در مدل رواناب ذوب برف بود. مقدار فاکتور درجه-روز محدود شده (a)‏ برای 


1 Senzeba et al. 

2 Nunchhani et al. 

3 Snowmelt Runoff Model 
4 Brubaker et al. 
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تمامی ماه‌های فصل ذوب در تمام سال‌های مورد مطالعه مقدار ۰/۲ برآورد شد. نتایج به‌دست آمده نشان داد 
که نسخه ترکیبی تابش عملکرد بهتری نسبت به نسخه درجه-روز متداول SRM‏ داشته است. وفاخواه و 
همکاران )۲۰۱٤(‏ مدل تابشی و درجه-روز SRM‏ را مقایسه کردند. آن‌ها در مدل تابشی SRM‏ از فاکتور 
درجه-روز محدود شده (,0) استفاده کردند و در این مطالعه. تابش خالص با استفاده از MODIS „oLa‏ 
برآورد شد. در نتایج به‌دست آمده مشخص شد ضریب فروکش و ضریب رواناب ذوب برف حساس‌ترین 
پارامترها هستند و مدل تابشی SRM‏ نتیجه بهتری با توجه به معیارهای دقت به‌دست آورد. استیلی و 
همکاران" (۲۰۱۷) محصولات برفی MODIS‏ برای مدل‌سازی رواناب ذوب برف را در حوضه ppl‏ رودخانه 
ریوگرانده مورد مطالعه قرار دادند. آن‌ها به کاربران توصیه می‌کنند هنگام شبیه‌سازی ذوب برف با مدل SRM‏ 
یکی از محصولات اصلاح شده با سایبان MODSCAG™‏ یا محصولات مخصوص به موقعیتشان را از 
مجموعه 0010۸1 Collection 6 NDSI‏ انتخاب و استفاده کنند. شی و همکاران" (۲۰۱۸) در مطالعه‌ای 
الگوریتم‌های بهینه‌سازی را برای کالیبراسیون پارامترهای متغیر در زمان مدل رواناب برف (SRM)‏ در شبیه- 
سازی جریان آب ناشی از ذوب برف برای حوضه رودخانه ماناسی بررسی کردند و نشان دادند برای شبیه- 
سازی یک فصل ذوب Gy‏ طول زیردوره SE‏ قابل‌توجهی بر عملکرد مدل دارد. نانج‌های و همکاران 
(۲۰۲۰) در پژوهشی تغییرات مکانی و زمانی پارامترهای برف در حوضه دیبانگ» آروناچال پرادش» هند را با 
استفاده از مدل رواناب ذوب برف توزیع شده فضایی ) (SDSRM‏ همراه با سنجش‌ازدور و تکنیک‌های GIS‏ 
را مورد تجزیه‌وتحلیل قرار دادند و نتیجه گرفتند که می‌توان آن را برای استفاده در حوضه‌های اندازه‌گیری 
نشده یا کمی اندازه‌گیری شده در منطقه هیمالیا توصیه کرد. زیمنس و همکارانش" (۲۰۲۱) تخلیه رواناب در 
حوضه بالادست رودخانه آتاباسکای در غرب کانادا را با استفاده از مدل رواناب ذوب برف SRM‏ و 
مشاهدات سنجش ازدور برف MODIS‏ را شبیه‌سازی کردند. نتایج نشان داد گنجاندن تغییرات مکانی و زمانی 
در فاکتور درجه-روز و جداسازی آن برای مناطق یخ‌زده و بدون توده برف برای بهبود شبیه‌سازی شرایط 
مختلف برف در طول چند سال مهم بوده است. همچنین آنها اظهار داشتند این مدل در شبیه‌سازی جریانات 
تاریخی به‌خوبی عمل کرده است. 

میریعقوب زاده و همکاران (۱۳۹۰) در حوضه آبخیز سد کرج نشان دادند که مدل SRM‏ با ضریب تبیین 
۷ و ۰/۹۶ به‌خوبی قادر به مدل کردن فرآیند رواناب حاصل از ذوب برف می‌باشد. در پژوهشی oS‏ 
ابراهیمی و همکاران (۱۳۹۵) با استفاده از مدل SRM‏ و تصاویر هشت روزه سنجنده MODIS‏ رواناب روزانه 


1 Steele et al. 

2 the MODIS Snow Covered-Area and Grain size retrieval 
3 Xie et al. 

4 Spatially Distributed Snowmelt Runoff Model 

5 Siemens et al. 


سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل SRM‏ و برآورد پارامتر فاکتور درجه-روز با استفاده از... ۵ 


حاصل از ذوب برف را برای حوضه آبریز نهاوند. مورد شبیه‌سازی و ارزیابی قرار دادند. بالاترین ضرایب 
همبستگی و تفاضل حجمی برای مدل در دوره واسنجی به‌ترتیب ۰/۷۵ و TW‏ و در دوره درستی‌سنجی 
پرابر ۰/۷۹ و ۲۰/۳۷ به‌دست آمد. تحلیل حساسیت پارامترهای ضریب درجه-روز. ضریب رواناب برف و 
باران مدل انجام گردید که OUS‏ دهنده حساسیت کمتر مدل به افزایش و کاهش ضریب رواناب باران در Ole‏ 
سایر پارامترها oy‏ اربابی سبزواری و نجاری ATAID‏ به پژوهشی مشابه در ارتفاعات حوضه آبخیز طالقان 
پرداختند. مدل» آبدهی روزانه را با مجذور ضریب همبستگی LEO‏ شبیه‌سازی کرد. نتایج این تحقیقات نشان 
می‌دهند که مدل شبیه‌ساز SRM‏ ابزاری بسیار مناسب برای برآورد میزان رواناب حاصل از ذوب برف می- 
باشد و از دقت بالایی برخوردار است و این مدل جهت مدیریت منابع آب کارآمد می‌باشد. در پژوهشی 
مشابه» رشیدی و همکاران (YAT)‏ با استفاده از مدل SRM‏ رواناب روزانه حاصل از ذوب برف را با 
استفاده از تصاویر سنجنده MODIS‏ در حوضه آبریز دربند سملقان محاسبه کردند. نتایج کاربرد مدل شبیه- 
سازی موفق و قابل قبولی را نشان داده است. نتایج برای تحلیل حساسیت مدل به پارامترهای ضریب رواناب 
برف» فاکتور درجه-روز و افت محیطی دما نشان داد که سه پارامتر یاد شده در بالا بر ضریب تبیین و درصد 
تفاضل حجمی مژثرند که البته مدل نسبت به تغییرات افت محیطی دما نسبت به دو پارامتر دیگر از حساسیت 
کمتری برخوردار است. در پژوهشی به‌منظور پایش زمانی و مکانی سطح پوشش برف. میزان ضریب تبیین 
(RO)‏ به‌دست آمده برای مدل‌های SRM‏ و HBV’‏ در مرحله واسنجی به‌ترتیب برابر ۰/۷۱ و ۰/۱۱ و مقدار 
۲ و ۰/۹ برای مرحله اعتبارسنجی به‌دست آمد (ارتیمانی و همکاران ۱۳۹۸). مدل SRM‏ به [Jo‏ استفاده 
از تصاویر ماهواره‌ای» نسبت به مدل HBV‏ عملکرد قابل قبول‌تری در شبیه‌سازی رواناب حاصل از ذوب 
برف از خود نشان داد. عالی‌نژاد و همکاران (۱۶۰۰) رواناب ناشی از ذوب برف را در حوضه‌ی قره‌سو با 
مدل SRM‏ شبیه‌سازی کردند. نتایج شبیه‌سازی بیانگر دقت بالای این مدل می‌باشد. به‌طوری‌که ضریب تبیین 
R?)‏ برای سال‌های آبی ۹۵-۹7 الی ٩-۹۷‏ به‌ترتیب Soles‏ ۰/۹۳ و ۰/٩‏ و درصد حطای حجمی (Dy)‏ آن نیز 
به‌ترتیب ۰/۳ و ۳/۳۳ به‌دست آمد. 

فقر آمار و اطلاعات در حوضه‌های آبریز یکی از چالش‌های جدی پیش روی مدل‌سازی رواناب ناشی از 
ذوب برف می‌باشد که کار تعیین و محاسبه پارامترهای ورودی به مدل را سخت و دشوار می‌سازد. یکی از 
پارامترهای ورودی مهم در مدل رواناب ذوب برف (SRM)‏ فاکتور درجه-روز (0) است. مقادیر فاکتور 
درجه-روز (0) در حوضه‌های آبریز به‌صورت متداول و کلاسیک توسط داده‌های اندازه‌گیری شده چگالی 
برف محاسبه و تعیین می‌شود و سپس به مدل SRM‏ وارد می‌گردد. این داده‌ها به دلیل توزیع مکانی و زمانی 
سطوح پوشیده از Gy‏ و کمبود ایستگاه‌های اندازه‌گیری در ارتفاعات مختلف دارای کیفیت پایین» کمبود و 
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فقدان در زمان و مکان‌های مختلف در حوضه‌های آبریز است. ازآنجاکه برف در مناطق صعب‌العبور تجمع 
می‌یابد. ایجاد ایستگاه‌های اندازه‌گیری و برف‌سنجی در ارتفاعات کاری سخت. زمان‌بر و پرهزینه است و با 
توجه به اینکه در مدل SRM‏ عمق ذوب برف روزانه رابطه مستقیم با پارامتر فاکتور درجه-روز و دما دارد و 
از سوی دیگر ad je‏ تابش خالص در محاسبات فاکتور درجه-روز (0) به روش متداول و کلاسیک آن در نظر 
گرفته نمی‌شود و در نتیجه تأثیر آن در عمق ذوب برف منظور نمی گردد» هدف از این مطالعه OF‏ است تا با 
استفاده از معادله‌ای جدید و برگرفته از معادلات گذشته تعیین عمق ذوب برف ارائه شده توسط کوستاس 
(۱۹۹۶) و بروبکر OAD‏ همراه با استفاده از تصاویر ماهواره‌ای تابش خالص و تکنیک‌ها سنجش‌ازدور به 
روشی ساده. کم‌هزینه دقیق‌تر و فیزیکی‌تر فاکتور درجه-روز (0) را محاسبه کرده و سپس polis‏ به‌دست 
آمده را جهت شبیه‌سازی رواناب ذوب برف برای مدل SRM‏ تعریف کرده و عملکرد روش جدید را در 
مقابل روش متداول و کلاسیک محاسبه فاکتور درجه-روز )0( در کنار لحاظ کردن polis‏ متفاوت پارامتر 
ضریب فروکش (O0‏ در شبیه‌سازی رواناب ذوب برف توسط مدل را ارزیابی و بررسی کنیم. در این مطالعه 
نتایج حاصل از شبیه‌سازی‌های مختلف رواناب ذوب برف توسط مدل SRM‏ با تعریف polis‏ فاکتور درجه- 
روز (A)‏ محاسبه شده با روش جدید و کلاسیک و تعریف polis‏ سالانه و فصلی برای پارامتر ضریب 
فروکش K)‏ مورد ارزیابی و مقایسه قرار گرفته است و بهترین روش محاسبه فاکتور درجه-روز (0) و 
بهترین شیوه شبیه‌سازی رواناب ذوب برف در حوضه آبریز مورد مطالعه توسط مدل SRM‏ تعیین گردیده 


۲- منطقه مورد مطالعه 


حوضه آبریز سد کارده با مساحتی حدود ۵۶6/۲۱۳ کیلومتر مربع در شمال شرق ایران و شمال استان 
خراسان رضوی و در شمال شهرستان مشهد قرار گرفته است. موقعیت جغرافیایی ol‏ بین Yo!‏ ?04 تا £o!‏ 
۳ طول جغرافیایی و ۳۸ ۳۹ تا oa"‏ ۵۷ ۳۷ عرض جغرافیایی قرار گرفته و این حوضه از زیر حوضه- 
های اصلی حوضه ppl‏ کشف‌رود به‌شمار می‌رود. کم ارتفاع‌ترین نقطه حوضه در محل سد کارده با ارتفاع 
۰ متر و بلندترین نقطه آن در ارتفاعات شمال خرکت با ارتفاع ۲۹۷۷ متر واقع شده است. 
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شکل ۱- موقعیت جغرافیایی منطقه و ایستگاه‌های مورد استفاده در حوضه آبریز کارده 


۳- مواد و روش‌ها 


مدل رواناب ذوب برف 4S (SRM)‏ توسط مارتینک 22 lavo m‏ توسعه داده سل یک روش مفهومی 


نیمه توزیعی درجه-روز cleo‏ برای شبیه‌سازی و پیش‌بینی رواناب روزانه حاصل از ذوب برف و بارندگی در 
es ts oes‏ ان مشک و OS SEM sale dol qa ates‏ ربب تشگنا کی 


درجه-روز و ضرایب رواناب برای برف و باران می‌باشند که می‌توانند شبیه‌سازی رواناب را اصلاح و تعدیل 
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کرده و بر عملکرد مدل تأثیر بگذارند (پاندی و همکاران؛ ۲۰۱۶). SRM‏ در گام زمانی روزانه و شبیه‌سازی 
جریان مبتنی بر تغیبرات دمای روزانه. بارش و پارامترهای مختلف به نمایندگی از شرایط فیزیکی حوضه 
عمل می XS‏ (الیاس و همکاران؛ ۵ فراهم کردن داده‌های عمده برای شبیه‌سازی و مدل‌سازی جریان در 
حوضه‌های کوهستانی و با توجه به نقش مهم رواناب ذوب برف و همچنین کمبود داده در این مکان‌ها یک 
چالش است (فیروزی و صادقیان. (XH VV‏ همچنین هرگونه تقویت در اساس فیزیکی روش درجه-روز استفاده 
از آن را در پیش‌بینی رواناب ذوب برف و ارزیابی اثرات تغییر اقلیم بهبود خواهد بخشید (کوستاس و رانگو" 
٤‏ 

مدل رواناب ذوب برف (SRM)‏ رواناب حاصل از ذوب برف و باران را به‌طور روزانه محاسبه کرده و با 
لحاظ کردن فروکش جریان» دبی خروجی از حوضه را برآورد می‌کند. مساحت و محدوده حوضه و ناحیه- 
cl‏ ارتفاعی بوسیله نقشه DEM.‏ منطقه و توسط نرم افزار ARCGIS‏ به‌دست آمد. منطقه به سه ناحیه 


ارتفاعی با فواصل ۰۰۰ متر تقسیم شد (جدول ۱). 


حدول ۱- مشخصات ناحیه‌های ارتفاعی تعیین شده در حوضه PI‏ کارده 


مساحت طبقات ار تفاعی ناحیه ارتفاعی 
(m) Km’ %‏ 
YMEE/AY VY‏ ۱۷۰۰۰ 
Yes ۳۳/۸۲ 0/40‏ 
VA/OY MURY‏ ۲۹۰۰-۲۰ 


] Panday et al. 
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4 Digital Elevation Model 


سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل SRM‏ و برآورد پارامتر فاکتور درجه-روز با استفاده از... 4 


59.425163 
E: h | E 
E: E 

راهنما 
طبقه بندی ارتقاع 
E‏ مرز-حوضه آبریز-کارده LL]‏ 4 & 
m 8‏ 1200-1700 ناحیه A‏ | 8 
m‏ 1700-2300 تاحیه B‏ `| 

= | MC تاحیه‎ 2300-2900 m = 
Re 0 3.75 7.5 15 Kilometers ۶ 
3 LE p xk P p x rk ] 3 

59.425163 


شکل ۲- نقشه طبقه‌بندی ارتفاعی حوضه آبریز کارده 


متوسط بارش روزانه از ایستگاه‌های هیدرومتری و باران‌سنجی موجود در منطقه برای سال‌های 
۲ ۱۱ و ۲۰۱۳ ۲۰۱۲ استفاده شد (جدول ۲) و توسط گرادیان بارش برای متوسط ارتفاع 
هییسومتریک هر aol‏ به‌دست آمد. دما متغیر مهمی در محاسبه عمق ذوب برف روزانه است. با توجه به 
پيشنهاد توسعه دهندگان مدل مبنی بر اينکه یک ایستگاه دمایی خوب (حتی اگر در خارج از حوضه قرار 
داشته باشد) به چندین ایستگاه کمتر قابل اطمینان ترجیح داده می‌شود. به این منظور از داده‌های دمای بيشینه 
و کمینه ایستگاه‌های سینوپتیک گلمکان و مشهد استفاده شد (جدول (Y‏ و توسط گرادیان دمایی برای متوسط 


ارتفاع هیپسومتریک هر یک از ناحیه‌های ارتفاعی استخراج شد. 
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حدول ۲- مشخصات ایستگاه‌های هیدرومتری و باران‌سنجی مورد استفاده در حوضه آبریز کارده برای سال‌های 
۱۲ و ۲۰۱۳ ۲۰۱۲ 


ارتفاع از سطح دریا طول حغرافیایی عرض جغرافیایی نوع ایستگاه نام ایستگاه 
(m)‏ 
Ya" °04 yq" ۱۳۳۲‏ ۳1 هیدرومتری کارده بالادست 
ory vy" ?o4 £ Y" ۱۳۷۹‏ باران‌سنجی سد کارده 
gy" 904 yq" IER‏ ۳1 باران‌سنجی آل 
9A °04 yo" 1۳‏ هیدرومتری و باران‌سنجی جنگ 
SY 2 904 vy" 1014‏ باران‌سنجی گوش بالا 
£A" 904 yy" ۱۸۳۰‏ ۳7 باران‌سنجی مارشک 
oe" °04 Yo" ۱۹۶۱‏ 7 باران‌سنجی بلغور 


حدول ۳- مشخصات ایستگاه‌های هواشناسی مورد استفاده در o>‏ ضه 3l‏ کارده برای سال‌های ۳ و 


1۲_11۳ 
ارتفاع از سطح دریا (m)‏ طول جغرافیایی عرض جغرافیایی نوع ایستگاه نام ایستگاه 
vé" 904 YA" 444‏ رت سینوپتیک مشهد 
ony va" 904 1۷" A‏ سینوپتیک گلمکان 


che‏ پوشش برف یکی از متغیرهای اصلی ورودی به مدل SRM‏ است که به‌صورت روزانه مورد نیاز 
مدل می‌باشد. برای استخراج سطح پوشیده از برف برای سال ۲۰۱۱۲۰۱۲ جهت واسنجی مدل و سال 
۳ برای درستی‌سنجی مدل از محصول پوشش برف روزانه سنجنده MODIS‏ از ماهواره Terra‏ 
به‌نام ۷10110۸1 استفاده گردید. در این محصولات برای استخراج سطح پوشش برف از الگوریتم snow‏ 
map‏ استفاده شده که در این الگوریتم از شاخص تفاضل نرمال شده پوشش برف ( ۷۳51 استفاده شده 
است. تصاویر به‌منظور تبدیل سیستم مختصات. عملیات برش تصاویر برای هریک از سه ناحیه ارتفاعی 
تعیین شده در حوضه آبریز کارده و محاسبه اطلاعات آماری به محیط نرم افزار ENVI‏ فراخوانده شد و 


] Normalised Difference Snow Index 
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شکل -Y‏ (الف) تغییرات سطح پوشش برف در سال ۱۲ و (ب) تغییرات سطح پوشش برف در سال 


۲۰۱۲.۳ در حوضه آبریز کارده 


ضریب رواناب برای کل حوضه ppl‏ محاسبه و در uo‏ مرحله واسنجی مقادیر بهینه آن به‌دست آمد. 
دمای بحرانی (Torr)‏ نیز Geb‏ توصیه توسعه‌دهندگان مدل, clea‏ ۲/۵ درجه سلسیوس برای ماه‌های سرد 
سال و دمای ۰/۷۵ درجه سلسیوس برای ماه‌های گرم سال در نظر گرفته شد. گرادیان دما (Y)‏ مقدار پیش 
فرض مدل و مقدار ۰/0۵ درجه سلسیوس به ازای هر ۱۰۰ متر ارتفاع و با توجه به روش انتخابی زمان تمرکز 
در مطالعات فیزیوگرافی که روش کرپیچ می‌باشده زمان تأخیر (L)‏ برابر VO‏ محاسبه گردید. 

فاکتور درجه-روز (0): فاکتور درجه-روز afem °c! d^]‏ تعداد درجه-روز Ted]‏ را به عمق ذوب 
برف روزانه [۷1]0 تبدیل می‌کند. فاکتور درجه-روز (0) در خارج از مدل محاسبه شده و سپس مقادیر آن 
برای مدل SRM‏ تعریف می گردد. فاکتور درجه-روز با در دسترس بودن مقادیر اندازه‌گیری شده چگالی 
برف» به شیوه کلاسیک و متداول از رابطه تجربی زیر به‌دست می‌آید. 

o=1.1% (\) 

که در این رابطه ,۵ چگالی برف و pu‏ چگالی آب می‌باشد. فاکتور ذوب به‌صورت مقدار عمق ذوب شده 

به تعداد درجه-روز بوده و به‌شکل رابطه (۲) بیان می‌گردد. در تلاشی بروبکر و همکاران GRE O81)‏ را در 

مدل ترکیب کردند. در این روش معادله ذوب برف به‌صورت رابطه (۳) بیان می‌شود. در این پژوهش با قرار 

دادن معادلات شماره (۲) و (۳) در یک دستگاه معادله جدیدی جهت برآورد فاکتور درجه-روز (0) به‌دست 
آمده است. معادله فاکتور درجه-روز (0) جدید به‌صورت رابطه شماره (۶) می‌باشد: 


M = a. Ta (Y) 
iu = mgo. Ra + ay. Ta (0) 


(£) 


NY‏ جغرافیا و مخاطرات محیطی شمارة اول 


mo. Rg 
= T, 


a + ar 


که در این معادله mg‏ ثابت فیزیکی تبدیل انرژی به حجم یا عمق آب است à, ۱/۰۲۱ em ۳/۷۷ day)‏ 
ضریب درجه-روز محدود شده است Ta .(cm/PC day)‏ شاخحص درجه Sole‏ است. Ra‏ تابش خالص 
روزانه است. به دلیل اینکه تغییرات فاکتور درجه-روز (0) به‌صورت ۱۵ تا ۳۰ روزه توصیه می‌شود برای 
محاسبه تابش خالص از محصول تاہبش خالص A‏ روزه NEO” (CERES -NETFLUX-E)‏ با قدرت تفکیک 
مکانی ۰/۲۵ درجه استفاده شد. جهت محاسبه GLU‏ خالص تصاویر در نرم افزار ENVI‏ فراخوانده شد و 
مقادیر آن استخراج گردید. با توجه به نتایج مطالعات مارتینک در سال ۱1۹۸۹ کوستاس و رانگو (1۹۹4) و 
بروبکر و همکاران (۱۹۹7) جهت برآورد مقدار ضریب درجه-روز محدود شده (Or)‏ در این مطالعه نیز مقدار 


۰/۲۵ برای این ضریب در محاسبات مورد استفاده قرار گرفت. 


جدول ۶- مقادیر پارامتر فاکتور درجه-روز (0) محاسبه شده با معادلات کلاسیک و جدید در حوضه کارده 


مقدار واقعی فرحه محدوده نرمال نماد نام متغیر 

"TET.‏ محاسبه فاکتور (0) با معادله کلاسیک 
cm ۰/۱-۱ ie‏ 

Jeca) : Di ۰/۲۵ — ۶‏ (« محاسبه فاکتور (0) با معادله جدید 


ضریب فروکش یک پارامتر بسیار مهم در مدل SRM‏ است. تحلیل داده‌های سری زمانی دبی اندازه‌گیری 
شده روش خوبی برای تعیین × می‌باشد Lb‏ توجه داشت که × ثابت نیست و با کاهش دبی مقدار OF‏ 
افزايش these‏ مدل SRM‏ خود K‏ را به‌عنوان ورودی نمی‌پذیرد و بایستی ضرایب ‏ و لبه مدل داده شود 
به این منظور مقادیر دبی جریان در هر روز Qua‏ در مقابل دبی روز قبل Qn‏ در قالب محور لگاریتمی رسم 


شد و با استفاده از نمودار بەدست آمده polio‏ ضرایب ‏ و y‏ محاسبه شد. 


حدول ۵- Cal pd polis‏ × و ¥ به‌دست آمده برای معادله فرو کش > E) SL‏ 


X Y 
ماه‌های سرد‎ * [AAA T 
ماه‌های گرم‎ */VA0Y (AV 
سالیانه‎ ۰/۱ YE 


1 Nasa Earth Observations 
2 The Clouds and the Earth's Radiant Energy System 


۱۳ برآورد پارامتر فاکتور درجه-روز با استفاده اتمه‎ s SRM سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل‎ ke 


مدل SRM‏ هیدروگراف رواناب اندازه‌گیری شده و محاسبه شده را به‌صورت گرافیکی نمایش می‌دهد. در 
عین حال برای ارزیابی عملی‌تر از چگونگی انجام شبیه‌سازی از دو معیار دقت» ضریب تبیین (RÊ)‏ و درصد 
تفاضل حجمی " (Dy)‏ استفاده می‌کند. مقادیر تعیین شده پارامترهای مورد نیاز مدل تا حد امکان براساس 
داده‌های موجود برآورد گردید و در صورت نبود اطلاعات و آمار کافی از مقادیر پیشنهادی توسعه‌دهندگان 
مدل و مطالعات مشابه استفاده و در مرحله واسنجی polis‏ بهینه آن‌ها تعیین شد. در این پژوهش مدل رواناب 
ذوب برف (SRM)‏ برای مقادیر به‌دست آمده فاکتور درجه-روز (0) از دو روش کلاسیک (معادله شماره ۱) 
و روش جدید doles)‏ شماره (E‏ در کنار مقادیر متفاوت ضریب فروکش جریان K)‏ به‌صورت یک مقدار 
ثابت برای plod‏ طول دوره شبیه‌سازی و استفاده از دو مقدار فصلی برای ماه‌های گرم و سرد اجرا شد و نتایج 
bol‏ برای Ged‏ بهترین روش شبیه‌سازی و عملکرد بهینه مدل توسط معیارهای ارزیابی دقت مدل مورد 


-٤‏ نتایج و بحث 


واسنجی مدل: مطابق آنچه شرح داده شد اطلاعات و داده‌های مورد نیاز برای اجرای مدل رواناب ذوب 
برف WInSRM‏ تعریف گردید و سپس مدل برای ماه‌هایی از سال که پوشش برف در سطح حوضه وجود 
dio cea eol‏ اب abo esae Tul Yo‏ اضرا مت ول colit‏ سای ۳۱۹۱۱۰۲۰۱۲ شام usas‏ 
مدل» متوسط بارش برف ۱/۷ سانتی‌متر در این سال است که نسبت به دهه قبل و بعد از آن بیشتر بوده و 
باعث عملکرد بهتر مدل در شبیه‌سازی رواناب ذوب برف می‌شود. شکل ۶. هیدروگراف‌های شبیه‌سازی شده 
و اندازه‌گیری شده GUL,‏ ذوب برف را در سال ۲۰۱۱-۲۰۱۲ که توسط مدل SRM‏ ترسیم شده را نشان 
می‌دهد. که برای ارزیابی و مقایسه بصری بهتر هیدروگراف‌های روش‌های مختلف در کنار یکدیگر آورده 
شده‌اند. براساس معیارهای ارزیابی دقت مدل (جدول CV‏ اجرای مدل رواناب ذوب برف با فاکتور درجه- 
روز جدید و استفاده از دو ضریب فروکش فصلی بهترین عملکرد را در شبیه‌سازی رواناب در طول دوره 
شبیه‌سازی داشت. در این اجراء ضریب تبیین (QD)‏ برابر با ۰/۷۲ و درصد اختلاف حجمی (Dy)‏ برابر با ۶/۱۷ 
بود. میزان حجم رواناب اندازه‌گیری شده ۷/۵۵۰ میلیون مترمکعب و متوسط رواناب اندازه‌گیری شده ۰/6۸ 
مترمکعب در انيه بود. شبیه‌سازی مدل در ماه‌های ژانویه, فوریه. مارس و آوریل از تطابق خوبی برخوردار 
بوده و جریان شبیه‌سازی شده به‌خوبی جریان مشاهده شده را دنبال می کند. مقدار حجم رواناب محاسبه شده 


۵ میلیون مترمکعب و متوسط روانات محاسبه شده ۰۸:۹ مترمکعب در asl‏ می‌باشد. مقادیر دبی و حجم 


] The coefficient of determination 
2 Volume difference 
3 Snowmelt Runoff Model for Windows 


M‏ جغرافیا و مخاطرات محیطی شمارة اول 


رواناب مشاهده شده و محاسبه شده بسیار به یکدیگر نزدیک هستند. پس از آن اجرای مدل با فاکتور درجه- 
روز کلاسیک و باز هم استفاده از دو ضریب فروکش فصلی در شبیه‌سازی موفق‌تر بوده و دارای ضریب Ox)‏ 
(R?)‏ برابر با ۰/0۲ و درصد اختلاف حجمی Dy)‏ برابر با VVÉA‏ بود. این روش دارای کم‌برآورد در مقدار 
OL >‏ شبیه‌سازی شده در ماه‌های نوامب دسامبر و فوریه و مارس بوده است. مدل SRM‏ ضعیف‌ترین 
عملکرد را در اجرای مدل با استفاده از فاکتور درجه-روز کلاسیک و یک ضریب فروکش سالانه برای کل 
دوره شبیه‌سازی را از خود نشان داد در این اجرا ضریب تبیین (RT)‏ و درصد اختلاف حجمی (Dy)‏ به‌ترتیب 
پرابر با ۰/6۳ و ۱۹/۲۲ بودند. در این روش مدل در ابتدای ماه نوامبر و ماه‌های پایانی دوره شبیه‌سازی دارای 


کم‌برآوردی می‌باشد و استفاده از این روش در شبیه‌سازی رواناب ذوب برف توصیه نمی‌شود. 


3 2 

3 1.5 

= 1 

1 0.5 

1 5 

5 ^y SN a ey y Vv 

۲ N N N N N M 

o e ER 5 SH SH‏ کي ت 
cial e RO d e‏ اندازه ی شده کی 


رواناب محاسبه شده با فاکتور Rp‏ و دو ضریب فروکش فصلی BEES‏ 
رواناب محاسبه شده با فاکتور کلاسیک و دو ضریب فروکش فصلی س  ..‏ 
رواناب محاسبه شده با فاکتور جدید و یک ضریب فروکش سالانه - - - - 
رواناب محاسبه شده با فاکتور »‌کلاسیک و یک ضریب فروکش سالانه ees‏ 


شکل -٤‏ رواناب اندازه‌گیری و شبیه‌سازی شده توسط مدل SRM‏ برای سال ۲۰۱۱_۲۰۱۲ در حوضه کارده 


جدول *- polis‏ معیارهای ارزیابی دقت مدل SRM‏ در اجراهای مختلف در سال ۲۰۱۱_۲۰۱۲ 


درصد اختلاف 
حجمی Dy‏ ضریب R? ond‏ روش‌های اجرای مدل 
NY &NW‏ رواناب محاسبه شده با فاکتور 0 جدید و دو ضریب فروکش فصلی 
AS 744۹‏ رواناب محاسبه شده با فاکتور 0 کلاسیک و دو ضریب فروکش فصلی 
VN‏ ۰/۱ رواناب محاسبه شده با فاکتور A‏ جدید و یک ضریب فروکش سالانه 
eY MUYY‏ رواناب محاسبه شده با فاکتور » کلاسیک و یک ضریب فروکش سالانه 


درستی سنجی مدل: شکل ۵ هیدر وگراف‌های tls)‏ اندازه گیری oJ‏ 3 شبیه‌سازی شده توسط مدل 
SRM‏ در حوضه آبریز کارده در سال ۲۰۱۳ ۲۰۱۲ که جهت صحت‌سنجی مدل اجرا شده را OUS‏ می‌دهد. 


همان‌طور که در (شکل Y‏ ب) مشاهده می‌شود در cnl‏ سال شاهد کاهش محسوس سطح پوشش برف و در 


سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل SRM‏ و برآورد پارامتر فاکتور درجه-روز با استفاده از o‏ 


نتیجه کاهش حجم آورد و نوسانات دبی مشاهده شده نسبت به دوره واسنجی هستیم که خود باعث کاهش 
lucida‏ ی رواب کوب ر ت فرط ال کرو اباب امین ماراق ارا کے 0 5یا 
eo d dise) ۴‏ ی( و ds‏ انادف disse d; D wee‏ با ووک 
محاسبه رواناب با فاکتور درجه-روز جدید و استفاده از یک ضریب فروکش سالانه به ترتیب برابر ۰/۵۱ و 
۸ بوده و دارای بهترین عملکرد در شبیه‌سازی رواناب توسط مدل SRM‏ از Ole‏ روش‌های موردبررسی 
در این تحقیق است و نتایج حاصل از دوره واسنجی را تأیید می‌کند و تطابق خوبی را برای دبی‌های اندازه- 
گیری شده و شبیه‌سازی شده در حوضه کارده OUS‏ می‌دهد. مدل در ماه‌های ابتدایی شبیه‌سازی, نوامبر تا 
نیمه دسامبر و ماه‌های پایانی شبیه‌سازی مارس و آوریل دارای کم‌برآورد بوده است و مقدار حجم رواناب 
دا کی 2 6۳۱ da gts caste dia‏ رالات اا کر de AA ad‏ مکی بر تاه ید 
باشد و مقادیر حجم رواناب محاسبه شده و متوسط رواناب محاسبه شده به‌ترتیب ۵/14۷ میلیون مترمکعب و 
۱ مترمکعب بر ثانیه است. پس از آن نیز اجرای مدل با فاکتور درجه-روز جدید و دو ضریب فروکش 
فصلی در شبیه‌سازی موفق‌تر بوده است. این اجرا دارای بیش‌برآورد در ماه‌های ژانویه. فوریه و آوریل "ue‏ 
و ریب a‏ ارو فرصت ات یی od E E ctus UD‏ سرا 
مدل با استفاده از روش محاسبه فاکتور درجه-روز کلاسیک و به‌کارگیری دو ضریب فروکش فصلی در طول 
روند شبیه‌سازی ضعیف‌ترین عملکرد مدل را از Ole‏ روش‌های به‌کار برده شده در این تحقیق به نمایش می- 
گذارد. در این اجرای مدل ضریب تبیین (Q2)‏ برابر با ۰/۳۳ و درصد اختلاف حجمی Dy)‏ برابر ۱۵/۳۲ می- 
باشد. در این شبیه‌سازی مدل دارای کم‌برآورد در ماه‌های نوامبر و دسامبر می‌باشد و نتایج حاصل از دوره 


زمان (روز) 0 
04/01/2013 03/01/2013 02/01/2013 01/01/2013 11/01/2012.:12/61/2002:0 
رواناب محاسبه شده با فاکتور »جدید و دو ضریب فروکش فصلی ج 
رواناب محاسبه شده با فاکتور ,کلاسیک و دو ضریب فروکش فصلی س . . — 
رواناب محاسبه شده با فاکتور جدید و یک ضریب فروکش سالانه - - - - 


شکل ۵- رواناب اندازه‌گیری و شبیه‌سازی شده توسط مدل SRM‏ برای سال ۲۰۱۱۲۰۱۲ در حوضه کارده 


bil ae n‏ و مخاطرات محیطی شماره ازل 


جدول polis - V‏ معیارهای ارزیابی دقت مدل SRM‏ در اجراهای مختلف در سال ۲۰۱۳ ۲۰۱۲ 


درصد اختلاف حجمی on wur‏ 


R? D,‏ روش‌های اجرای مدل 

rw‏ ۷ رواناب محاسبه شده با فاکتور 0 جدید و دو ضریب فروکش فصلی 
jtm ۱۹/۳۳‏ رواناب محاسبه شده با فاکتور 0 کلاسیک و دو ضریب فروکش فصلی 
1۳۸ \0/+ رواناب محاسبه شده با فاکتور 0 جدید و یک ضریب فروکش سالانه 
71۲ ۳۷ رواناب محاسبه شده با فاکتور 0 کلاسیک و یک ضریب فروکش سالانه 


کشور ایران در منطقه گرم و خشک قرار دارد و میزان حجم توده‌های برفی در حوضه‌های آبریز این گونه 
مناطق مانند حوضه آبریز کارده محدود می‌باشد و این یکی از دلایل کاهش دقت شبیه‌سازی توسط مدل 
است. میزان متوسط بارش برف در حوضه آبریز کارده در سال ۲۰۱۲_۲۰۱۳ برابر با ۳۰/۰۵ سانتی‌متر بوده و 
از دلایل کاهش دقت مدل در دوره درستی‌سنجی نسبت به دوره واسنجی می‌توان به کاهش ۲۷/۰۲ درصدی 
بارش Gy‏ و کاهش آورد سالانه اشاره کرد. همچنین خطای انسانی در اندازه‌گیری رواناب را از ps‏ عمده 
ol‏ دانست. چنانکه قربانی‌زاده و همکاران (۱۳۸۹) نیز در مطالعه خود که بر روی رواناب حاصل از ذوب 
برف با کمک مدل SRM‏ در حوضه آبریز کارون انجام دادند و دی و همکاران" (AA)‏ در شبیه‌سازی ذوب 
برف در پاکستان به عدم تطابق دبی‌های حداکثر مشاهده شده و شبیه‌سازی شده اشاره کردند و بیان کردند که 
این ضعف ناشی از اشتبامات در داده‌های مشاهداتی و بی‌کیفیتی آن‌ها می‌باشد. نتایج به‌دست آمده از شبیه- 
سازی رواناب حاصل از ذوب برف در این مطالعه با نتایج حاصل از تحقیقات میریعقوب‌زاده و همکاران 
(YA)‏ اربابی سبزواری و نجاری (Y)‏ و آرتیمانی و همکاران (۱۳۹۸) مطابقت دارد. همچنین مطالعه 
قاسمی و همکاران (۱۳۹۲) که نوسانات رواناب حاصل از ذوب برف تحت SU‏ پدیده تغییر اقلیم در دهه- 
های آینده را موردبررسی قرار دادند و در واسنجی مدل ضریب تبیین و تفاضل حجمی را به‌ترتیب ۰/۷۱ و 
۵ درصد به‌دست آوردند نتایج به‌دست آمده را تأیید می‌کند. آن‌ها نتایج حاصل از شبیه‌سازی واسنجی 
مدل را قابل قبول ارزیابی کردند. نتایج مطالعه لی و ویلیامز" (۲۰۰۸) بر روی رواناب حاصل از ذوب در 
حوضه بدون آمار تاریم در چین که مقدار ضریب تبیین را در دوره واسنجی ۰/۶ و برای دوره درستی‌سنجی 
۸ و ۰/۵۱ به‌دست آورده‌اند نیز با نتایج تحقیق حاضر هم‌خوانی دارد. نتایج حاصل از محاسبه فاکتور 
درجه-روز با استفاده از روش جدید (معادله شماره (E‏ و روش محاسبه کلاسیک (معادله شماره (Y‏ و مقایسه 
نتایج آنها نیز با مطالعات مارتینک A8)‏ کوستاس و رانگو OIE)‏ و بروبکر و همکاران )144( وفاخواه 
و همکاران QE)‏ زیمنس و همکاران (۲۰۲۱). دلاور و همکاران (۱۳۹۰). ابراهیمی و همکاران (۱۳۹۵). 


1 Dey et al. 
2 Li and Williams 
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رشیدی و همکاران (YT)‏ مطابقت دارد. نتایج حاصل از اجرای مدل با روش‌های مختلف برآورد فاکتور 
درجه-روز در دوره‌های واسنجی و درستی‌سنجی SE‏ از برتری روش جدید برآورد فاکتور درجه-روز 
(معادله شماره (E‏ در شبیه‌سازی رواناب ناشی از ذوب برف نسبت به روش کلاسیک محاسبه آن (معادله 
شماره (V‏ است. نتایج به‌دست آمده OUS‏ می‌دهد با استفاده از روش جدید محاسبه فاکتور درجه-روز و 
محصول تابش خالص هشت روزه NEO (CERES-NETFLUX-E)‏ می‌توان با روشی ساده‌تر و کم‌هزینه به 
مقادیر قابل قبولی از فاکتور درجه‌-روز )0( دست پیدا کرد. همچنین نتایج تحقیق حاضر OUS‏ می‌دهد 
ضریب فروکش جریان (K)‏ پارامتری حساس در شبیه‌سازی رواناب ذوب برف توسط مدل 58۷ است و 
تأثیر به‌سزایی در تطابق جریان اندازه‌گیری شده و شبیه‌سازی شده دارد و در نتیجه تعیین دقیق Ol‏ می‌تواند 
نتایج بهتری با توجه به سنجه‌های کارایی مدل به‌دست آورد. نتایج حاصل از تحقیقات وفاخواه و همکاران 
(۲۰۱۶) نیز این مطلب را تأیید می‌کند. نتایج حاصل OUS‏ می‌دهد با توجه به اینکه مقدار ضریب فروکش 
جریان متغیر بوده. در دوره‌هایی که نوسانات جریان و آورد رودخانه 5 236 است مدل SRM‏ می‌تواند با استفاده 
از دو مقدار ضریب فروکش hed‏ در شبیه‌سازی نقاط اوج موفق‌تر عمل XS‏ و در دوره‌های که OU m‏ 
یکنواخت بوده و نوسانات آن کم است. استفاده از یک ضریب فروکش سالانه عملکرد بهتری دارد. 


۵- نتیجه گیری 


با توجه به نتایج حاصل از این مطالعه و اجرای مدل با روش‌های مختلف در دوره‌های واسنجی و 
درستی‌سنجی. بهترین روش مدل‌سازی جریان در دوره واسنجی» روش محاسبه فاکتور درجه-روز جدید 
(معادله شماره (E‏ و کاربرد دو مقدار ضریب فروکش جریان فصلی با ضریب تبیین ۰/۷۲ و درصد اختلاف 
حجمی ۶/۱۷ بود. در دوره درستی‌سنجی نیز روش شبیه‌سازی رواناب با فاکتور درجه-روز جدید (معادله 
شماره (E‏ و یک ضریب فروکش سالانه بهترین نتایج را از نظر معیارهای ارزیابی دقت مدل ارائه داد. ضریب 
تبیین در OF‏ ۰/۵۱ و درصد اختلاف حجمی 1/۳۸ بود. ازاین‌رو می‌توان نتیجه گرفت محاسبه مقادیر فاکتور 
درجه-روز (0) با استفاده از معادله جدید ابداع شده در این مقاله (معادله شماره (E‏ و تعریف مقادیر آن برای 
مدل SRM‏ به مدل کمک می‌کند تا عملکرد بهتری نسبت به شیوه کلاسیک محاسبه فاکتور (0) با استفاده از 
داده‌های چگالی برف اندازه‌گیری شده در حوضه ppl‏ (معادله شماره ۱) ارائه دهد. ترکیب مؤلفه تابش در 
معادله محاسبه فاکتور (0) نیز نتایج مطالعات گذشته مانند کوستاس و رانگو (۱۹۹۶) و بروبکر و همکاران 
)1441( در ترکیب مولفه تابش در محاسبه عمق ذوب برف در مدل SRM‏ را تأیید می‌کند. محصول تابش 
خالص استفاده شده در اين مطالعه برای محاسبه فاکتور )0( به روش جدید NEO (CERES-NETFLUX-‏ 


E)‏ نیز از دقت بسیار خوبی برخوردار بوده است و می‌تواند به نحو قابل قبولی برای شیبه‌سازی رواناب 


\A‏ جغرافیا و مخاطرات محیطی شمارة اول 


حاصل از ذوب برف در حوضه‌هایی که با کمبود داده مواجه هستند مورد استفاده قرار گیرد. همچنین نتایج 
معیارهای ارزیابی دقت ارائه شده نشان می‌دهد ضریب فروکش جریان نقشی اساسی در تطابق جریان‌های 
اندازه‌گیری شده و شبیه‌سازی شده Lal‏ می‌کند و همانطور که در مطالعات گذشته مانند وفاخواه و همکاران 
(۲۰۱۶) آمده پارامتری حساس و موثر در شبیه‌سازی رواناب ذوب برف توسط مدل SRM‏ است. نتایج 
ع ail als QUSE tag cul‏ شرب مر اند VN Loy oe‏ فان راخ U5 4 SRM div‏ 35$ 
و تعیین دقیق آن در بهبود عملکرد مدل کمک شایانی می‌کند. ضریب فروکش جریان مقداری ثابت نیست و 
با کاهش دبی مقدار OF‏ افزایش می‌یابد و بالعکس. شبیه‌سازی رواناب ناشی از ذوب برف در حوضه آبریز 
کارده نشان می‌دهد برف قسمت قابل توجهی از آورد رودخانه در فصل زمستان و بهار را تشکیل می‌دهد. 
ذوب برف از اواخر ماه فوریه در حوضه افزایش می‌یابد و می‌تواند در صورت همراه شدن با رواناب حاصل 
از بارش باران منجر به وقوع سیل به حصوص در اواخر ماه مارس و ماه آوریل گردد که نیازمند برنامه‌ریزی و 
مدیریت کارآمد در جلوگیری از خسارات ناشی از سیلاب‌ها و استفاده بهینه از این منابع آبی می‌باشد. با 
توجه به نتایج این پژوهش مطالعه بر روی رواناب ناشی از ذوب برف در حوضه‌های آبریز امری ضروری 
است و می‌تواند در مدیریت بحران و برنامه‌ریزی برای مقابله با مخاطرات و خسارت‌های حاصل از سیلاب- 
ها les SETS REN‏ وماق و ااا قفا عار نک هه وی ills‏ وان 
کرده و از بروز خسارات جبران‌ناپذیر جلوگیری به عمل آورد. اين‌گونه مطالعات بر روی زمان‌بندی و ميزان 
رواناب حاصل از ذوب برف به مدیران در برنامه‌ریزی و مدیریت کارآمد برای تأمین و توزیع آب مورد ناز 
بخش کشاورزی در فصل‌های گرم سال. جهت جلوگیری از مخاطرات ناشی از کمبود آب بەخصوص در 
سال‌هایی که کم‌بارش بوده یا به علت GES‏ زودهنگام هوا ذوب زودهنگام توده‌های برف اتفاق می‌افتد. 
کمک شایانی می‌نماید. با توجه به اینکه کشور ایران در منطقه خشک قرار گرفته مطالعات مستمر این‌گونه 


منابع برای مدیریت مؤثر و کارآمد و جلوگیری از هدررفت آب ضروری می‌باشد. 
این مقاله از رساله دکتری با عنوان ارزیابی اثرات تغییر اقلیم بر توده‌های برف و رواناب سطحی (مطالعه 
موردی: حوضه آبریز کارده) در دانشگاه فردوسی مشهد استخراج شده است. 
کتاینامه 
— ارتیمانی» محمد مهد ES‏ زینی‌وند. em‏ طهماسبی‌پور: ناصر؛ AYIA‏ ارزیابی کارآیی مدل SRM‏ و HBV‏ در -dui‏ 


سازی رواناب ناشی از ذوب برف در حوزه آبخیز بوجین. تحقیقات مناب عآب ایران» ۵( ۲۶۱-۲۲۸. 
https://www.sid.ir/paper/100221/fa‏ 
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